ROTATIONAL MODULATION OF THE RADIO EMISSION 
FROM THE M9 DWARF TVLM 513-46546: BROADBAND 
COHERENT EMISSION AT THE SUBSTELLAR BOUNDARY? 

G. Hallinan 1 , A. Antonova 2 , J.G. Doyle 2 , S. Bourke 1 , W. F. Brisken 3 and A. Golden 1 

ABSTRACT 

The Very Large Array was used to observe the ultracool rapidly rotating 
M9 dwarf TVLM 513-46546 simultaneously at 4.88 GHz and 8.44 GHz. The 
radio emission was determined to be persistent, variable and periodic at both 
frequencies with a period of ~ 2 hours. This periodicity is in excellent agreement 
with the estimated period of rotation of the dwarf based on its v sin i of ~ 60 km 
s _1 . This rotational modulation places strong constraints on the source size of 
the radio emitting region and hence the brightness temperature of the associated 
emission. We find the resulting high brightness temperature, together with the 
inherent directivity of the rotationally modulated component of the emission, 
difficult to reconcile with incoherent gyrosynchrotron radiation. We conclude that 
a more likely source is coherent, electron cyclotron maser emission from the low 
density regions above the magnetic poles. This model requires the magnetic field 
of TVLM 513-46546 to take the form of a large-scale, stable, dipole or multipole 
with surface field strengths up to at least 3kG. We discuss a mechanism by which 
broadband, persistent electron cyclotron maser emission can be sustained in the 
low density regions of the magnetospheres of ultracool dwarfs. A second non- 
varying, unpolarized component of the emission may be due to depolarization 
of the coherent electron cyclotron maser emission or alternatively, incoherent 
gyrosynchrotron or synchrotron radiation from a population of electrons trapped 
in the large-scale magnetic field. 

Subject headings: radio continuum: stars - stars: low-mass, brown dwarfs - 
stars: activity - stars: magnetic fields - stars: rotation - radiation mechanisms: 
non-thermal 
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1. INTRODUCTION 

Magnetic activity in cool main sequence stars is usually diagnosed through the detection 
of coronal X-ray and radio emission, and via the emission from spectral lines such as chro- 
mospheric Ha. In the Sun's case, it is accepted that such magnetic activity is powered by 
magnetic fields generated by a so-called afl dynamo at the base of the convective layer where 
differential rotation is strongest (Parker 1975). This dynamo mechanism is also thought 
to operate in a wide spectral range of cool stars from F-type to early M-type dwarfs and 
is particularly successful in accounting for the well-established rotation-activity relationship 
observed over this spectral range (Pallavicini et al. 1981; Noyes et al. 1984). 

For low mass stars with spectral type > M3, the nature of magnetic activity is less 
well understood. Here lies the boundary at which low mass stars become fully convective, 
no longer possessing the radiative layer necessary to sustain the aVL dynamo. However, M 
dwarf flare stars, very often with spectral type later than M3, are associated with intense 
levels of magnetic activity and are also characterized by large magnetic filling factors, with 
kG magnetic fields covering a large proportion of the stellar disk (Saar 1994; Johns-Krull & 
Valenti 1996). Indeed, the fraction of objects exhibiting emission associated with magnetic 
activity, particularly high levels of quiescent Ha emission, is seen to increase steadily for 
spectral types > M3 reaching 70% by about spectral type M8 (Gizis et al. 2000; West 
et al. 2004). Recently, magnetic mapping of the stellar surface of an M4 dwarf has directly 
confirmed the presence of large-scale, axisymmetric, poloidal fields, including a strong dipolar 
component, despite the absence of surface differential rotation (Donati et al. 2006). Clearly 
an extremely efficient alternative dynamo operates in low mass, fully convective stars. 

Beyond spectral type ~ M8, the realm of ultracool stars and brown dwarfs, a sharp drop 
occurs in the levels of quiescent Ha and X-ray emission observed, along with a complete 
breakdown of the activity-rotation relationship, with many fast rotators (v sin i > 15 km 
s _1 ) showing no evidence of quiescent Ha or X-ray emission (Neuhaser et al. 1999; Gizis et 
al. 2000; Mohanty & Basri 2003; Fleming et al. 2003; West et al. 2004). This reduction 
in activity may be due to the reduced ionization levels in the atmospheres of these cooler 
dwarfs (Mohanty et al. 2002). Magnetic fields become decoupled from the surrounding 
atmosphere and are no longer able to generate the necessary magnetic stress required for the 
dissipation of heat to sustain a chromosphere and corona. However, this does not preclude 
the efficient operation of a dynamo or the generation of large-scale kG magnetic fields on 
ultracool dwarfs. In fact, despite the reduction in the levels of quiescent emission, there have 
been detections of Ha and X-ray flares from a number of ultracool dwarfs (Reid et al. 1999; 
Gizis et al. 2000; Rutledge et al. 2000; Liebert et al. 2003; Fuhrmeister & Schmitt 2004; 
Stelzer 2004; Rockenfeller et al. 2006a) including a T dwarf (Burgasser et al. 2000) which 
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point to the presence of strong magnetic fields. Rockenfeller et al. (2006b) have recently used 
multiband photometric monitoring to confirm a periodicity of 3.65 ±0.1 hours in observations 
of the M9V dwarf 2M1707+64. Moreover, these authors associate this periodicity with the 
rotational period of the dwarf and rule out the presence of dust clouds as a source of this 
variability. Instead, they attribute this periodicity to the presence of magnetically induced 
cool spots on the surface of the star, further evidence for the presence of strong surface 
magnetic fields on ultracool dwarfs. 

Based on the sharp drop in quiescent Ha and X-ray emission and the empirical relation- 
ships established between the observed X-ray and radio emission for a wide spectral range of 
low mass, main sequence stars (Giidel & Benz 1993; Benz & Giidel 1994), it was expected 
that any radio emission from ultracool stars and brown dwarfs would be below the sensitivity 
of the current generation of radio telescopes. Nonetheless, Berger et al. (2001) unexpectedly 
detected persistent, highly variable radio emission from the M9 field brown dwarf LP 944-20 
at levels which, based on the lower limit for the quiescent X-ray flux established by the 
Chandra observation of Rutledge et al. (2000), violated the Gudel-Benz relations by more 
than 4 orders of magnitude. Follow up VLA observations by Berger (2002, hereafter B02) 
yielded three further candidates. All detected sources with measured rotation rates were 
fast rotators, with confirmed v sin i ranging up to 60 km s -1 , despite quiescent levels of Ha 
and X-ray emission being absent, greatly suppressed, or not yet observed. B02 postulated 
that the rot at ion- activity relationship may hold for radio emission in ultracool stars and 
brown dwarfs despite breaking down in other wavelength regimes. Burgasser & Putman 
(2005) also conducted a survey of southern late M and L dwarfs using the Australia Tele- 
scope Compact Array (ATCA) and confirmed two detected sources. One of these sources, 
DENIS 1048-3956, was seen to emit an extremely bright (~ 20 mJy), narrow bandwidth, 
100% circularly polarized flare, with the inferred brightness temperature > 10 l3 K requiring 
a coherent emission mechanism, most likely an electron cyclotron maser which would require 
the presence of kG magnetic fields. Berger et al. (2005, hereafter B05) have also confirmed a 
periodicity of ~ 3 hours in the radio emission from the L3.5 candidate brown dwarf 2MASS 
J00361617+1821104 (hereafter 2M 0036+18) in observations separated by three years, citing 
three possible sources of the periodicity, (1) orbital motion of a companion (2) rotation of 
2M 0036+18 or (3) weak periodic flaring. Finally, a large scale survey of late M, L and 
T dwarfs has yielded three more confirmed radio detections and highlights the trend that 
Lrad/Lhoi increases with later spectral type, contrary to what is observed at Ha and X-ray 
wavelengths (Berger 2006). 

Observations thus far have established that persistent radio emission at the bottom of 
the main sequence is a reality, despite the sharp drop in activity observed in X-ray and 
Ha emission and confirming the probable existence of large-scale, stable magnetic fields on 
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ultracool stars and brown dwarfs. In this paper we report on VLA observations of a previ- 
ously detected M9 dwarf, TVLM 513-46546 (hereafter TVLM 513), undertaken to further 
investigate the mechanism responsible for this anomalous radio emission. These observations 
provide insights into the nature, topology and strength of the associated magnetic fields and 
more fundamentally, the dynamo mechanism responsible for their generation. 



2. TVLM 513 



TVLM 513 is a young disk M9 dwarf located at a distance of 10.5 pc with a surface 
temperature of 2200 K and a bolometric magnitude of logLb o i/L ~ —3.65 (Tinney et 
al. 1993, 1995; Leggett et al. 2001). The absence of lithium in the spectrum of TVLM 513 
requires its mass to be > 0.6 M Q (Reid et al. 2002). We further loosely constrain its mass to 
be < 0.8 M for an age r < 10 Gyr, based on correlation of its determined properties with 
the evolutionary models put forward by Baraffe et al (1998) and Chabrier et al. (2000). This 
places TVLM 513 right at the substellar boundary, being either a high mass brown dwarf or 
an older, very low mass star. It is one of the most rapidly rotating ultracool dwarfs thus far 
detected with a measured v sin i of ~ 60 km s _1 (Basri 2001). Despite this rapid rotation it 
displays only weak Ha: emission, with equivalent width measurements of 1.7-3.5 A (Martin et 
al. 1994; Reid et al. 2001; Mohanty & Basri 2003), and there are no reported X-ray detections 
in the literature. Although Ha activity is suppressed, it has been strongly detected at radio 
frequencies. B02 used the VLA to detect persistent, variable emission at 8.46 GHz which 
included a strongly right circularly polarized (~ 65%) flare which reached a flux intensity of 
~ 1100/iJy. More recently, Osten et al. (2006) conducted a multi-frequency VLA observation 
of TVLM 513 at 8.4 GHz, 4.8 GHz and 1.4 GHz, using a strategy that involved time-sharing 
a single 10 hour observation between the various frequency bands. TVLM 513 was detected 
at each frequency band with only marginal confirmation of variability and no detection of 
flares or strong circular polarization, possibly due to the absence of continuous coverage at 
any one frequency band. 



3. RADIO OBSERVATIONS 

The radio observations of TVLM513 were conducted with the NRAO Very Large Array 
(VLA) 1 , for a duration of 4.74 hours on 2005 January 13 and 4.74 hours on 2005 January 



1 Thc VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 
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14, simultaneously in C (4.88 GHz) and X (8.44 GHz) band. This was achieved by splitting 
the VLA into two subarrays, with 14 and 11 antennas used for C and X band respectively 
on January 13, and 13 and 12 antennas used for C and X band respectively on January 14. 
3C286 was used for flux calibration and 1513+236 was used for phase calibration with a time 
on source in a single scan of 12 minutes before moving to the phase calibrator for 2 minutes 
for both bands. Reduction of the data was carried out using the AIPS software package. 
After applying the standard procedures of calibrating and editing the data, the background 
sources were removed from the visibility data. This was achieved by using the task IMAGR 
to CLEAN the region around each source, and then the task UVSUB to subtract the resulting 
models of the background sources from the map. The task UVFIX was used to shift the 
tangent point coordinates of the target source to coincide with the phase centre. Light curves 
were generated by plotting the real part of the complex visibilities at the position of TVLM 
513 as a function of time. 



4. PROPERTIES OF THE RADIO EMISSION 

TVLM 513 was detected as a variable source at both 4.88 GHz and 8.44 GHz (Figure 
1). The average flux received was 405 ± 18 /iJy at 4.88 GHz and 396 ± 16 yuJy at 8.44 GHz, 
with the corresponding spectral index of a = 0.0 ±0.2 consistent with observations of earlier 
spectral type M dwarfs (Giidel & Benz 1996). There was significant variability around these 
average flux levels, as highlighted in Figure 1, which plots the total intensity (Stokes I) radio 
light curves received at each frequency for both observations with a time resolution of 15 
minutes. In order to quantify this variability we assumed a constant source and determined 
the reduced chi-square values, xl-, f° r Stokes I radio light curves derived from the data at 
each frequency with a range of time resolutions from 5-60 minutes. The resulting xl values 
range from 1.2 - 4.0 for the 4.88 GHz light curves and 1.4 - 4.8 for the 8.44 GHz light curves. 
The higher xl values are generally derived from the light curves with lower time resolution, 
due to reduced error in the data points as a result of lower rms noise in the corresponding 
longer integration images. Based on the maximum Xr values attained, the 4.88 GHz data 
agree at the 4.63% level with a constant source while the 8.44 GHz data agree at the 2.9% 
level with a constant source. The upper limit (3a) on the net degree of circular polarization 
is r c < 13% at 4.88 GHz and r c < 12% at 8.44 GHz. However, in higher time resolution 
light curves there are significant (> 3cr) detections of both left and right circularly polarized 
emission at both frequencies. Also, as mentioned above, B02 reported onasi 65% right 
circularly polarized flare at 8.44 GHz. As we will discuss further below, we postulate that 
there may be components of both left and right circularly polarized emission from TVLM 
513 which may add destructively to give a much reduced net circular polarization. 
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The radio emission from TVLM 513 is also periodic at both frequencies with a period 
of ~ 2 hours. This is particularly clear in Figure 2 which correlates the Stokes I radio 
flux received at 4.88 GHz on January 13th and 14th 2005. To investigate this periodicity we 
performed a Lomb-Scargle periodogram analysis (Lomb 1976; Scargle 1982; Press & Rybicki 
1989) of the 4.88 GHz and 8.44 GHz Stokes I radio light curves for a wide range of time 
resolutions. The significance of the peaks derived from this analysis was dependent on the 
choice of time resolution, with the strength of the periodic signal reduced for more coarsely 
sampled data. With this in mind, the highest time resolution available of 10 seconds is used 
for both the 4.88 GHz (Figure 3) and 8.44 GHz (Figure 4) periodograms. Two clear peaks 
are present in both periodograms at frequencies of ~ 0.5 hr -1 and ~ 1 hr -1 with the peak at 
0.5 hr" 1 corresponding to the fundamental 2 hour period. The secondary harmonic peak at 
1 hr -1 is accounted for by the presence of two main peaks in the emission per 2 hour period 
(Figure 1). The manifestation of both peaks as a group of multiple peaks (spectral leakage) is 
due to the large gap (« 17 hours) in the observation. We used a program provided by Harry 
Lehto (private communication) which used a 1-dimensional CLEAN algorithm (Roberts et 
al 1987) to iteratively remove the spectral window function from the raw data and confirm 
the presence of both peaks. To further assert that this periodicity is not an artifact of the 
observation, we conducted a Lomb-Scargle periodogram analysis of the 4.88 GHz and 8.44 
GHz Stokes I light curves of a field source detected in the observation of TVLM 513. A time 
resolution of 10 seconds was used and the field source data have the same spectral window 
function as the TVLM 513 data. No significant peaks were present in the periodogram of 
the field source data. 

In order to determine the significance of the detected periodicity in the radio emis- 
sion from TVLM 513 we conducted Monte Carlo simulations producing Lomb-Scargle peri- 
odograms of 60,000 randomized versions of the data at both frequencies. For all randomized 
datasets the same spectral window function as the original data was also maintained. The 
resulting false alarm probabilities are shown in Figure 3 and Figure 4. Both peaks are 
present with > 99% probability in the 4.88 GHz data and > 99.9% probability in the 8.44 
GHz data. We also include the periodogram of the data from the single observation at 8.44 
GHz on January 14 2005 (Figure 5). Once again, the peaks at ~ 0.5 hr -1 and ~ 1 hr -1 arc 
present with > 99.9% probability. It is also notable that the spectral leakage observed in 
the periodograms in Figure 3 and Figure 4 is absent in this periodogram which is derived 
from a single observation and does not contain a large gap in the data. Interestingly, there 
are also further significant peaks present in the periodogram at the frequencies of 2 hr -1 , 3 
hr -1 and 5 hr -1 , harmonics of the fundamental 0.5 hr -1 frequency. This may be evidence 
of further structure in the periodic emission from TVLM 513, in particular the periodic 
light curve may be non-sinusoidal in nature. Similar harmonic peaks were not detected in 
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the periodograms of the other individual datasets and further observations are required to 
investigate this possibility. 

The data were epoch folded to generate a high signal to noise light curve of the periodic 
radio emission from TVLM 513. We initially folded the data for a range of possible periods 
at or near two hours and analysed the Xr values of the resulting light curves in an attempt 
to detect a peak in variance. A clear peak was observed at approximately 1.96 hours and 
this value was used as the putative period for epoch folding the data. Phase bins of 6, 7 and 
8 minutes were used with the resulting light curves plotting the total intensity (Stokes I) 
and circular polarization (Stokes V) of the data shown in Figure 6 (4.88 GHz) and Figure 7 
(8.44 GHz). For each figure two phases of the epoch folded light curve are shown for clarity. 

At both frequencies the emission is shown to vary in a non-sinusoidal manner with 
periodic sharp increases in intensity characterized by high levels of either left or right circular 
polarization. At 8.44 GHz there is a short periodic right circularly polarized burst where 
total intensity levels exceed 1000 /zJy and circular polarization levels reach 40%. It is worth 
noting that a similar burst was observed in a previous 1.83 hour, 8.44 GHz observation of 
TVLM 513 (B02). At 4.88 GHz a broad peak in emission observed at phase xs 0.3 is seen to 
consist of two adjacent but separate peaks, the first showing strong left circular polarization 
and the second showing right circular polarization. This amounts to a periodic reversal in 
the polarization of the emission. To investigate the variation in circular polarization with 
frequency we also correlate high time resolution light curves of the Stokes V, 4.88 GHz and 
8.44 GHz data, shown in Figure 8. The large periodic peak in left circular polarization 
observed in the 4.88 GHz data has no corresponding counterpart in the 8.44 GHz data, 
possibly indicating a polarization reversal from low to high frequencies. Another possibility 
is that different regions producing emission with opposite polarities are contributing to the 
emission and different regions dominate at different frequencies. The presence of various 
regions producing either left or right circularly polarized emission would also account for the 
periodic reversals in the polarization of the emission at 4.88 GHz. 

We can therefore assert that the radio emission from TVLM 513 is characterised by 
periodic bursts of both left and right circularly polarized emission. It is noteworthy, how- 
ever, that the flux intensity does not drop below ~ 200/iJy at either frequency and circular 
polarization is absent for much of the epoch folded light curve. This indicates that the 
radio emission from TVLM 513 consists of two components, a highly circularly polarized 
periodically varying component and a non- varying, unpolarized component. 

Periodicity was not detected previously by B02, as the corresponding 1.83 hour ob- 
servation was shorter than the estimated period of ~ 2 hours. Similarly, the observation 
conducted by Osten et al. (2006) involved time-sharing between 3 frequency bands with 



- 8- 



the coarseness of the resulting time-series probably inhibiting the detection of the periodic 
signal. 

5. SOURCE OF THE PERIODICITY 

The most plausible explanation for the periodicity of the radio emission from TVLM 
513 is its period of rotation. The maximum period of rotation can be constrained using the 
measured v sin i of ~ 60 km s _1 and an estimated radius, the latter being a strong function 
of both mass and age. Based on a T e s of 2200 K, the inferred age is greater than 400 Myr, 
which together with a mass estimate of 0.6 M — 0.8 M yields a radius of ~ 0.1 ± 0.01 R 
(Chabrier et al. 2000), which is in agreement with the estimated radius for TVLM 513 of 
Dahn et al. (2002). The resulting maximum period of rotation is 2.0 ± 0.2 hours. This is 
in excellent agreement with the observed periodicity of the radio emission from TVLM 513 
and infers a high inclination angle i > 65°, i.e., the equatorial plane is very close to our line 
of sight. Considering the striking similarities in the properties of the radio emission from 
TVLM 513 and the L3.5 dwarf 2M J0036+18 (B02; B05), which include periodicity, a high 
degree of variability and the periodic presence of high levels of circular polarization, it seems 
probable that the same emission mechanism is responsible in both cases. Although, B05 
cited 3 possible sources for the 3 hour periodicity in the radio emission from 2M 0036+18, 
rotation of the dwarf would seem the most likely candidate in light of the probable rotational 
modulation of the radio emission from TVLM 513. Based on the v sin i of 15 ± 5 km s _1 
measured for 2M J0036+18 by Schweitzer et al. (2001), this would require a low inclination 
angle i 24 ± 8° of the axis of rotation, i.e. almost pole-on. However, Zapatero Osorio 
et al. (2006) have recently reported on a revised measurement of 36 ± 2.7 km s _1 , in close 
agreement with the equatorial velocity derived from a 3 hour period of rotation. This would 
require a much larger inclination angle between the axis of rotation and our line of sight 
> 66°, similar to what is found for TVLM 513. 

Any viable model for the radio emission from these two ultracool dwarfs must account 
for the rotational modulation, which coupled with the high degree of variability, place strong 
constraints on the length scale of the emitting regions and therefore the brightness temper- 
ature of the associated radio emission. 
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6. RADIO EMISSION MECHANISMS 

To date, the radio emission from TVLM 513, and other detected ultracool stars and 
brown dwarfs, has been attributed to incoherent gyrosynchrotron emission from a non- 
thermal population of mildly relativistic electrons (B02; Burgasser & Putman 2005; B05; 
Osten et al. 2006), with the exception of the highly polarized flare detected from DENIS 
1048-3956 (Burgasser & Putman 2005). Such a population of electrons is thought to follow 
a power law distribution, n(E) oc E~ 5 , with 5 « 2 — 4 typically invoked for low mass stars 
(Gudel 2002). The resulting radio emission is broadly peaked perpendicular to the field 
(f] u oc gin ^~°- 43+0 - 65<5 , Dulk & Marsh 1982), which, if applicable to the periodic component 
of the radio emission from TVLM 513, requires occultation of emitting active regions by the 
stellar disc to account for the variability and periodicity. Indeed, the high degree of variabil- 
ity of the radio emission, characterized by three-fold variations in flux over short timescales 
< 0.5 hours (Figure 6 & Figure 7), requires the distribution of the emitting active regions 
to be non-uniform and confines their length scale, L, to be much less than the size of the 
stellar disc. 

The brightness temperature associated with an incoherent radiation process is limited 
to ~ 10 12 K by inverse Compton cooling (Kellerman & Pauliny-Toth 1969; Readhead 1994). 
However, this limiting value applies to synchrotron emission from ultra-relativistic electrons 
in very weak magnetic fields <C 1 G, as observed from extragalactic sources, whereas here 
we are interested in gyrosynchrotron emission from mildly relativistic electrons and the 
corresponding higher strength stellar magnetic fields. In this case, the emission is emitted at a 
frequency, v, where 10—100 and the electron cyclotron frequency v c ps 2.8 x 10 6 5 

Hz. Dulk & Marsh (1982) have shown that the effective temperature of gyrosynchrotron 
emission, T e g, from an isotropic distribution of electrons with a power-law index 2 < 5 < 7 
and a low energy cut-off of 10 keV is 

T cS » 2.2 x io9lO- - 315 (sin^)- a36 - ao65 (-) a50+ao85<5 K (1) 

and is limited to a few times 10 9 K. The brightness temperature for all incoherent emission 
is also limited by T B < T eS (Dulk 1985), consequently T B is thus limited to a few times 10 9 
K for gyrosynchrotron emission. 

The brightness temperature of the radio emission from TVLM 513 is given by 



T B = 2 x 10 9 (f u /mJy)(v/GHzy\d/pc)\L/R Jup y 2 K 



(2) 
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Assuming a radius R ps 0.1 R ~ Rj up infers a T B ~ 3.7 x 10 9 (L/i? Jup )~ 2 , based on the 
average flux received at 4.88 GHz of 405 /zJy. Even assuming a length scale L of the order 
of the size of the stellar disk would yield a brightness temperature ~lx 10 9 K. Indeed, this 
source was also detected at 1.4 GHz by Osten et al. (2006) with a brightness temperature 
of Tb =~ 2.9 x 10 10 (L/i?j up )~ 2 . If gyrosynchrotron emission is applicable, such a high 
brightness temperature requires a source size of the order of the stellar disk or greater. 

Contrary to this, the characteristics of the periodic light curve require highly directive 
emission from compact regions much less than the size of the stellar disk with a corresponding 
much higher brightness temperature. In particular, the highly circularly polarized bursts at 
both frequencies (Figure 7 & Figure 8) have very narrow duty cycles (< 0.15 phase) occurring 
over timescales much less than that required for an emitting region to rotate in an out of 
view, thus confirming the directivity of the emission. Moreover, the periodic reversal of the 
sense of polarization of the emission at 4.88 GHz confirms the presence of more than one such 
active region. It seems unlikely, therefore, that the high brightness temperature and high 
directivity of the rotationally modulated component of the radio emission from TVLM 513 is 
compatible with isotropic gyrosynchrotron emission from a large extended corona, although 
it cannot be ruled out as a possible source of the non- varying, unpolarized component. 

Other incoherent emission processes such as synchrotron emission from a population of 
ultra relativistic electrons can be highly beamed perpendicular to the local magnetic field, 
accounting for the variability and rotational modulation of the radio emission, and can reach 
higher brightness temperatures than gyrosynchrotron emission. Indeed, synchrotron radia- 
tion has been previously proposed as a source of rotationally modulated stellar radio emission 
(Lim et al. 1992, 1994). However, the detection of significant levels of circular polarization 
is incompatible with synchrotron emission. Alternatively, a coherent emission process may 
be applicable with the possibility of much higher brightness temperatures being invoked. 
One such coherent process is plasma radiation, which is generated at the plasma frequency, 

1/2 

Up ps 9000n e Hz where n e is the plasma electron density. However, plasma radiation is more 
dominant at lower frequencies (< 1 GHz) due to increasing free-free absorption with increas- 
ing v (and hence n e ) (Giidel 2002). It has been speculated that the degree of absorption 
may be relaxed for a higher temperature, dense, coronal plasma enabling the generation of 
plasma radiation at higher frequencies (White & Franciosini 1995; Osten & Bastian 2006). 
However, such conditions are at odds with those expected for the increasingly cooler and 
more neutral atmospheres of ultracool dwarfs as highlighted by the reduced levels of X-ray 
and Ha emission. 

Another possible coherent emission process is electron cyclotron maser emission which, 
as discussed in §7, can account for the properties of the rotationally modulated emission 
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from both TVLM 513 and 2M 0036+18, particularly if the magnetosphere takes the form of 
a large-scale, stable, dipole or multipole, with kG field strengths at the surface, such as that 
confirmed for earlier type M dwarfs by Donati et al. (2006). 

7. ELECTRON CYCLOTRON MASER EMISSION 

The electron cyclotron maser (hereafter referred to as maser) is a particle-wave plasma 
instability caused by the resonance between the gyrating electrons in an externally gener- 
ated magnetic field and the electric field of electromagnetic waves at frequencies near the 
electron cyclotron frequency and perhaps its low harmonics. It has been proposed as a 
possible generation mechanism for a wide range of astrophysical radio emission including 
planetary radiation from all of the magnetized planets, solar microwave bursts as well as 
certain emissions from dMe stars, T Tauri stars, RS CVn binaries, Algol-like binaries and 
magnetic chemically peculiar stars. 

To date, masering has been overlooked as a possible source of the radio emission from 
TVLM 513 and 2M 0036+18 due to the broadband, persistent nature of the detected emission 
from both candidates (B05, Osten et al. 2006), whereas maser emission is expected to be 
confined to short timescales and narrow bands around the fundamental or second harmonic 
of the electron cyclotron frequency of the emitting source region. However, the emission 
is only narrow-banded if the masering is confined to a region with very little variation in 
the local magnetic field strength and hence, the electron cyclotron frequency. This may be 
expected if the anisotropic electron distribution responsible for the maser instability solely 
takes the form of a loss-cone, a model applied to terrestrial auroral kilometric radiation by 
Wu & Lee (1979). In this model, electrons trapped in a magnetic field mirror at the surface, 
developing a loss-cone distribution due to the precipitation of electrons with small pitch 
angles. This results in an excess in the perpendicular component of the electron velocity 
distribution perpendicular to the magnetic field 5u± such that 5f /5v± > 0, providing the 
free energy to power the maser. This model was further developed and suggested as a source 
of certain solar and stellar bursts by Melrose & Dulk (1982), who proposed that mirroring 
electrons trapped in coronal loops may form a loss-cone at the loop footpoints. In the case 
of maser emission powered by a loss-cone, the instability is often confined to regions close 
to the stellar or planetary surface limiting the range in magnetic field strength and hence 
the bandwidth of the emission. It is also subject to saturation of the wave growth as the 
radiative process removes the loss-cone anisotropy from the electron distribution. 

However, work by Pritchett (1984a,b), Pritchett & Strangeway (1985) and Winglee & 
Pritchett (1986) highlighted an alternative mechanism by which the electron distribution 
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can also achieve the necessary anisotropy required for maser emission. Electrons accelerated 
along magnetic field lines into regions of higher field strength conserve magnetic momen- 
tum by evolving adiabatically to higher pitch angles, generating the necessary excess in the 
perpendicular component of the electron velocity distribution perpendicular to the magnetic 
field, 5v±. This "shell distribution" of electrons was subsequently determined to be more 
dominant in the generation of terrestrial auroral kilometric radiation than the loss-cone dis- 
tribution (Ergun et al. 2000, and references therein). Interestingly, maser emission from 
a shell distribution of electrons is not confined to regions close to the stellar or planetary 
surface, but can occur from source regions at any height where conditions are suitable for 
the operation of the maser. Also, the saturation conditions which limit the growth and 
continuous generation of maser emission from a loss-cone distribution do not apply to the 
shell maser, which can sustain steady-state, high brightness temperature emission as long 
as electron acceleration in the source region is maintained (Pritchett et al. 1999). The con- 
ditions in the source region must be such that (a) the magnetic field strength is relatively 
high and the plasma electron density resulting in an electron cyclotron frequency v c much 
greater than the plasma frequency u p , where v c m 2.8 x 10 6 B Hz and v v « 9000ny 2 Hz, and 
(b) electrons are continuously accelerated in the source region and can adiabatically evolve 
to form the shell distribution required to power the maser. If these conditions were present 
over a large range of heights above the surface of ultracool dwarfs, with a corresponding wide 
range of magnetic field strengths, maser emission could viably be continuously generated at 
a wide range of frequencies resulting in broadband, persistent radio emission. 

This model, if applicable to ultracool dwarfs, is a significant departure from the model 
generally applied to cool stars, which attributes the bulk of broadband, persistent radio 
emission from such stars to incoherent, gyrosynchrotron radiation. In fact, the coherent 
radio emission observed from all of the magnetized planets in our solar system (Zarka 1998, 
and references therein), attributed to the electron cyclotron maser instability, is of more 
relevance to this discussion. In particular, we highlight the coherent radio emission from the 
quasi-dipolar magnetosphere of Jupiter, which originates in the low density, high magnetic 
latitude regions and from the magnetic flux tube of the volcanic moon Io. This radio emis- 
sion is emitted in a number of components at a range of heights above the planetary surface, 
with frequencies ranging up to a maximum of approximately 40 MHz which corresponds 
to the electron cyclotron frequency at the regions of strongest magnetic field strength at 
Jupiter's surface (approx 14 G). The low cutoff frequency corresponds to the electron cy- 
clotron frequency at the height where v c /v p falls below a critical value, typically between 
~ 2.5 and 10, thereby quenching the instability (Zarka et al. 2001). The emission propagates, 
for the most part, in the supraluminous X mode, resulting in 100% circularly or elliptically 
polarized emission highly beamed perpendicular to the local magnetic field, that can reach 



-13- 



brightness temperatures up to 10 20 K. This X mode emission has an inherent right-handed 
circular polarization relative to the local magnetic field, resulting in the detection of right 
circularly polarized emission from Jupiter's northern hemisphere and left circularly emission 
from Jupiter's southern hemisphere. Certain components of this radio emission, such as the 
decametric emission (~ 10 to 40 MHz) from the high magnetic latitudes, referred to as non- 
Io DAM, can be broadband (Az/ ~ u), non- varying for timescales of the order of minutes, 
and show strong rotational modulation. 

By comparison, we now consider this maser emission process in the context of the 
radio emission properties and magnetospheric conditions on ultracool dwarfs such as TVLM 
513 and 2M 0036+18. It has been established that the magnetic fields are large-scale and 
stable due to the presence of periodicity in radio observations of 2M 0036+18 separated 
by two years (B05). If this large scale, stable magnetosphere takes the form of a dipole 
or multipole, with kG fields present at the surface there may be a large proportion of the 
magnetosphere where source conditions are suitable for the generation of maser emission. In 
particular, considering the atmospheres of dwarf stars become increasingly neutral beyond 
about spectral type M7 (Mohanty et al. 2002) it is reasonable to assume that the average 
plasma electron density also drops significantly. This assumption is supported by the sharp 
drop in detectable levels of Ha and X-ray emission for late M and L dwarfs. Moreover, if 
the magnetosphere takes the form of a dipole or multipole, the plasma density would not 
be homogeneous throughout the magnetosphere. For example, in the case of a dipole, any 
plasma injected into the magnetosphere would remain trapped mirroring in the closed field 
lines and possibly form "radiation belts" analogous to Earth's Van Allen belts. However, 
in regions over the poles, plasma would rapidly escape along open field lines leading to the 
formation of density cavities in the plasma. In this scenario, the plasma density would be 
much higher in the radiation belts than in the "coronal holes" above the magnetic poles. 
The rapid rotation of the dwarf might also further aid the trapping of particles at lower 
latitudes. It is feasible that the plasma density in these density cavities is such that v v <C v c 
for a wide range of heights above the stellar surface. Assuming emission at the fundamental 
of the electron cyclotron frequency, the upper cutoff frequency of emission would probably 
correspond to the electron cyclotron frequency at the stellar surface while the lower cutoff 
frequency of emission would correspond to the electron cyclotron frequency at the height 
above the surface where v c jv v is insufficient and the instability quenches. In the case of 
TVLM 513 and 2M 0036+18, with simultaneous detections at frequencies of 4.9 GHz and 
8.5 GHz, this would require a range of magnetic field strength in the density cavity ranging 
from 1.7 - 3 kG with a plasma electron density C 3 x 10 11 cm~ 3 . As discussed in §1 the 
probable detection of an electron cyclotron maser flare at a frequency of 8.5 GHz from the 
M8.5 dwarf DENIS 1048-3956 (Burgasser k Putman 2005) would indicate that kG field 
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strengths may be present on ultracool dwarfs. A dipolar magnetic field can be modelled by 



where B s and R s are the surface magnetic field strength and radius respectively. Considering 
that the radius of an evolved ultracool dwarf is approximately the same as that of Jupiter, 
the length scale of the emitting region producing simultaneous emission at 4.9 GHz and 8.5 
GHz on TVLM 513 and 2M J0036+18 would be smaller than that producing broadband 
decametric emission in the frequency range ~ 10 — 40 MHz from Jupiter's polar regions, 
albeit with much stronger kG fields in the source region. An electron density of <C 3 x 10 11 
cm" 3 in the low density regions above the magnetic poles of an ultracool dwarf is extremely 
likely, and may be many orders of magnitude less than this value, considering coronal electron 
density estimates for more massive, hotter, X-ray emitting stars range from 10 s cm" 3 — 10 13 
cm" 3 (Giidel 2004). If by some mechanism, a component of the population of trapped 
mirroring electrons was continuously accelerated into these low density plasma cavities due 
to the presence of magnetic field-aligned electric fields, these regions could be the source of 
persistent, broadband, coherent emission. The generation of such quasi-stable, long-lived, 
electric fields would have to be inherent to the nature of the large-scale, stable magnetosphere. 
One possible source of electric potential drops may be the formation of electrostatic double 
layers between the populations of trapped ions and electrons. Ions and electrons mirroring at 
different locations could possibly form a strong magnetic field-aligned potential drop, which 
may be a source of the continuous particle acceleration necessary for steady-state maser 
emission. 

The resulting radiation is expected to be highly beamed perpendicular to the magnetic 
field, explaining the rotational modulation of the emission from both TVLM 513 and 2M 
J0036+18. For 2M J0036+18 the emission shows strong net left circularly polarization 
indicating that the emission from one hemisphere is dominant. Indeed, Figures 3 & 4 in B05 
would seem to indicate that this circular polarization reaches almost 100% once per period 
of emission. In the case of TVLM 513, efficient masering from both hemispheres is required 
to account for the properties of the radio emission, with detections of both left and right 
circularly polarized emission. This is unsurprising considering the rotational axis is almost 
perpendicular to our line of sight. The presence of two strong peaks per period of emission 
separated by ~ 0.5 phase also strongly argues for the presence of a dipolar component to the 
magnetosphere. The broad peaks in emission occur when the magnetic axis of the dipolar 
field is perpendicular to our line of sight, which occurs twice per rotational period. However, 
the fact that there is no net circular polarization for much of its rotational phase requires 
a rather contrived geometrical coincidence if the emission is generated solely in the density 
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cavities associated with a dipolar field. The presence of a multipolar component to the 
magnetic field, with the corresponding associated density cavities also producing coherent 
radio emission, may account for this low net polarization. It is also possible that mode 
conversion or depolarization of the X mode emission occurs in the density cavity. Mode 
conversion of terrestrial auroral kilometric radiation from X mode to R mode in the emitting 
density cavity is proposed by Ergun et al. (2000) as an efficient method by which the maser 
emission can escape without reabsorption at higher harmonics of the emission frequency. 
Similarly, mode conversion from X mode to O mode has also been considered as a possible 
source of the elliptical polarization of the Jovian decametric emission (Zarka 1998, and 
references therein). 

Alternatively, the non- varying unpolarized component may be incoherent gyrosynchrotron 
or synchrotron radiation. Indeed, if the magnetic field takes the form of a large-scale, stable 
dipole or multipole, the population of electrons mirroring in closed field lines may con- 
tribute to the radio emission. As discussed above, this population of electrons would have 
an anisotropic distribution which, for gyrosynchrotron emission, can result in higher inten- 
sity emission than that produced by an isotropic population (Fleishman & Melnikov 2003). 
It is also possible that such a population of electrons may be continuously accelerated by 
radial diffusion to extremely high ultrarelativistic energies conducive to the generation of 
synchrotron emission. Such emission is generated at decimetric wavelengths from the radia- 
tion belts of Jupiter and may be a contributing component to the emission from both TVLM 
513 and 2M J0036+18. 



8. DISCUSSION AND CONCLUSIONS 

We report on VLA observations showing rotational modulation of the radio emission 
from the M9 dwarf TVLM 513, which places strong constraints on the source size of the 
emitting region, and hence brightness temperature of the emission. The resulting high 
brightness temperatures along with the implicit directivity of the rotationally modulated 
component of the radio emission are difficult to reconcile with incoherent gyrosynchrotron 
radiation from a population of mildly relativistic electrons. We discuss a new model for the 
radio emission from ultracool dwarfs based on the assumption that the magnetic field on 
these objects takes the form of a large scale dipole or multipole powered with kG strength 
fields at the surface. In regions over the magnetic poles, plasma can escape along open 
field lines forming density cavities where the electron cyclotron frequency is much greater 
than the plasma frequency. If a component of the population of trapped electrons was 
continuously accelerated into such density cavities by magnetic field-aligned electric fields, 
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possibly due to the presence of electrostatic double layers, such density cavities may be the 
source of steady-state, broadband electron cyclotron maser emission. We also discuss the 
possible presence of a non-varying, unpolarized component of the radio emission, which may 
be due to depolarization or mode conversion of the maser emission in the density cavity. 
Alternatively, this component of the radio emission may be incoherent gyrosynchrotron or 
synchrotron radiation. In particular, trapped particles mirroring on closed field lines in 
the large scale magnetic field can form high density radiation belts where electrons can be 
accelerated to ultrarelativistic energies emitting synchrotron radiation. 

Coherent maser emission can account for the presence of high levels of circular polariza- 
tion and the high brightness temperatures of the radio emission from ultracool dwarfs. The 
periodicity and high degree of variability from both TVLM 513 and 2M J0036+18 can be 
attributed to the inherent directivity of the radio emission coupled with the rapid rotation 
of the dwarf. We differentiate between this broadband, persistent maser emission and the 
extremely bright 100% circularly polarized flare observed from the M8.5 dwarf DENIS 1048- 
3956 (Burgasser & Putman 2005), typical of flares observed from earlier type M dwarfs. 
Such flares are almost always narrow-banded, 100% circularly polarized and can display 
structure over very short millisecond timescales. We speculate that both the broadband, 
quiescent emission and narrowband, flare emission are generated in the same source regions 
by the same mechanism, the electron cyclotron maser, albeit with fundamental differences 
in the nature by which the emission is generated. For example, the smooth broadband emis- 
sion may be associated with constant particle acceleration due to quasi-stable electric fields 
across the source region, whereas the narrowband flaring may be due to more impulsive 
acceleration events. This is analogous, once again, to Jupiter's decametric emission, which 
is characterized by smooth, broadband bursts with Av ~ v (referred to as L bursts), and 
more impulsive, narrowbanded bursts with Av <C v, which can reach much higher brightness 
temperatures (referred to as S bursts) (Zarka 1998). 

The configuration of the magnetosphere may also be a contributing factor to the ubiqui- 
tous rapid rotation of ultracool dwarfs. It is unlikely that such a large-scale stable magnetic 
field would undergo reconnection events in the increasingly neutral atmospheric plasma of 
ultracool dwarfs, minimizing momentum loss due to dissipation of magnetic energy, and thus 
leading to much longer spin down times than observed on earlier type dwarfs. This model 
for the radio emission from ultracool dwarfs can also account for the fact that L ra d/£boi 
is maintained, and even seen to increase for late M and L type dwarfs, despite the sharp 
drop in magnetic activity detected at other frequencies. The lower plasma density with later 
spectral type which inhibits the generation of Ha and X-ray emission, leads to conditions 
where more efficient generation of radio emission by the electron cyclotron maser instability 
can occur. 
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The mechanism by which Ha and X-ray flares are generated on ultracool dwarfs is 
still unclear. Chabrier & Kiiker (2006) postulate that the magnetic stresses necessary for 
such flares are generated in the higher temperature conductive layers of the star and rise to 
the surface, dissipating magnetic energy in the upper atmospheric layers. Such a model is 
supported by the observations of Fuhrmeister & Schmitt (2004) who report on a Ha flare 
from the M8.5 dwarf DENIS 1048-3956, the source of the large coherent radio flare discussed 
in §1. They interpret broadening on the blue side of the spectral line during the flare as 
possible evidence that the flare was associated with a rising cloud of emission ejected from 
the dwarf. Such ejections may also be the source of the trapped plasma mirroring in the 
large-scale magnetic field responsible for the generation of the radio emission. 

The generation of coherent emission in density cavities associated with large scale dipolar 
or multipolar fields should be considered as a possible source of both quiescent and flaring 
radio emission in a wide range of astrophysical objects, including T Tauri stars (Smith et 
al. 2003), RS CVn binaries (Osten et al. 2004), Algol- type binaries (Mutel et al. 1998) and 
chemically peculiar stars (Trigilio et al. 2000). Moreover, the absence of high levels of circular 
polarization and narrowband structure does not preclude the possibility that such emission 
may be coherent. 

The recent confirmation of a strong, large-scale axisymmetric field on a rapidly rotating 
M4 dwarf, with a strong dipolar component (Donati et al. 2006), highlights the possibility 
that electron cyclotron maser emission from polar density cavities may be a viable source of 
radio emission from dMe stars. Indeed, Bingham et al. (2001) have suggested maser emission 
from a shell distribution of electrons as a possible source of the 100% circularly polarized 
emission from UV Ceti. Considering the ubiquity of highly polarized, narrowband, coherent 
bursts from such objects, it is quite possible that broadband, slowly varying maser emission 
from the same source regions responsible for the flares may also contribute to the quiescent 
component of the emission. Indeed, such objects have much longer periods of rotation 
than later spectral type dwarfs such as TVLM 513 and 2M J0036+18 resulting in much 
slower rotational modulation of the highly beamed emission, and hence lower variability. 
An important factor in the determination of the emission mechanism in dMe stars has 
been the direct measurement of the extent of stellar radio coronae, and hence constraint 
of the brightness temperature of the associated emission, afforded by very long baseline 
interferometry (VLBI). Some VLBI observations of single dMe stars have revealed large 
extended radio coronae up to a few stellar radii in diameter, with brightness temperatures 
and degrees of circular polarization in line with incoherent emission such as gyrosynchrotron 
or synchrotron radiation (Alef et al. 1997; Benz et al. 1998; Pestalozzi et al. 2000). However, 
other observations of single dMe stars have revealed unresolved flaring and quiescent emission 
with brightness temperatures, Tb > 10 10 K, and very high degrees of circular polarization 
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(Benz & Alef 1991; Benz et al. 1995). In these cases a coherent emission mechanism such 
as electron cyclotron maser is the more plausible candidate. Further VLBI observations 
may clarify the degree to which both coherent and incoherent mechanisms contribute to the 
radio emission from dMe stars, and whether the prevalence of an individual mechanism is 
frequency dependent, activity dependent or related to the magnetospheric conditions on each 
individual star. 

Further high sensitivity observations of ultracool dwarfs previously detected at radio 
frequencies should be undertaken to investigate the presence of periodicity and to confirm if 
this periodicity is indeed due to the rotation of the dwarf. Direct detection of the magnetic 
fields on these objects may also be feasible. Assuming a gyrosynchrotron source would 
imply a magnetic field strengths of the order of a few hundred gauss to lkG, however, 
the coherent process, requires considerably higher field strengths of several kilo-gauss which 
could potentially be detectable using the FeH band around l//m (Reiners & Basri 2006). 
Although this technique is still in its infancy, the above authors show that that are sufficient 
FeH lines available even for the faster rotators whose observation and modelling could have 
an important bearing on the question of which radio process is at work in these ultracool 
dwarfs. 

Confirmation of the generation of broadband coherent emission from ultracool dwarfs 
would have immense implications for our understanding of stellar magnetic activity and 
the nature of the dynamo mechanism responsible for the generation of the magnetic fields 
driving this activity. It would also highlight the possibility that a similar mechanism for the 
generation of coherent radio emission from large-scale dipolar or multipolar fields may be a 
operating in a wide range of astrophysical objects, ranging from planets to brown dwarfs and 
cool stars. Comparison with more exotic phenomena such as coherent pulsar radio emission 
is also compelling. 
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Fig. 1. — Light curves for the Stokes I radio flux received from TVLM 513 at 4.88 GHz (top) 
and 8.44 GHz (bottom) on 2005 January 13 and 2005 January 14. The time resolution used 
is 15 minutes for all light curves. 
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Fig. 2. — Light curves for the Stokes I radio flux received from TVLM 513 at 4.88 GHz on 
2005 January 13 (solid line) and 2005 January 14 (dashed- dotted line). The two light curves 
have been aligned to highlight the close correlation and periodicity. 
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Fig. 3. — Lomb-Scargle periodogram of the 4.88 GHz Stokes I radio emission from TVLM 
513 with a time resolution of 10 seconds. Clear peaks are present at 0.5 hr -1 and 1 hr -1 , 
corresponding to the two peaks in emission per 2 hour period. The spectral leakage present 
for each peak is due to the spectral window function (shown on top) which is a result of 
the 17 hour gap between the two observations. The small peak in the spectral window 
function at ~ 4 hr -1 is due to the phase calibration every 15 minutes. The dashed-dotted 
lines represent the false alarm probabilities of the significance of the peaks based on 60,000 
Monte-Carlo simulations of randomised versions of the light curve with the same spectral 
window function as the original data. The solid lines on the right hand side of the graph 
represent the false alarm probabilities of the significance of the peaks as calculated by the 
Lomb-Scargle algorithm. 
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Fig. 4. — Same as Figure 3 for the 8.44 GHz Stokes I radio emission from TVLM 513 with 
a time resolution of 10 seconds. Once again, two clear peaks at 0.5 hr _1 and 1 hr _1 are 
present. 
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Fig. 5. — Lomb-Scargle periodogram of the 8.44 GHz Stokes I radio emission from the single 
observation of TVLM 513 on 2005 January 14, with a time resolution of 10 seconds. There 
is no spectral leakage due to the absence of a large gap in the data. The single peak in 
the spectral window function at ~ 4 hr" 1 is due to the phase calibration every 15 minutes. 
Once again, the two clear peaks at 0.5 hr -1 and 1 hr" 1 are present. Significant peaks at the 
frequencies of 2 hr -1 , 3 hr -1 and 5 hr -1 , harmonics of the fundamental 0.5 hr -1 frequency, 
are also present. This highlights the possible presence of further structure in the periodic 
radio light curve of TVLM 513, in particular the periodic light curve may be non-sinusoidal 
in nature. 
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Fig. 6. — Epoch folded light curves derived from the 4.88 GHz Stokes I (top) and Stokes V 
(bottom) radio data from TVLM 513 with phase bins of 6, 7 & 8 minutes. Two periods of 
the same epoch folded light curve are shown for clarity. Positive values for Stokes V indicate 
right circular polarization while negative values indicate left circular polarization. The main 
peak at phase 0.3 is resolved to consist of two adjacent peaks, the first one showing strong 
left circular polarization with the second showing right circular polarization. These highly 
polarized emissions add destructively over the course of one period of emission resulting in 
much lower net circular polarization. 
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Fig. 7. — Epoch folded light curves derived from the 8.44 GHz Stokes I (top) and Stokes 
V (bottom) radio data from TVLM 513 with phase bins of 6, 7 & 8 minutes. Two periods 
of the same epoch folded light curve are shown for clarity. Positive values for Stokes V 
indicate right circular polarization while negative values indicate left circular polarization. 
The presence of periods of right circular polarization is clear while evidence for left circular 
polarization is marginal at best. Of particular interest is the short peak of right circularly 
polarized emission where Stokes I flux levels exceed 1000 /iJy. The presence of such short 
duration periodic peaks implies that the emission is inherently directive or beamed in nature. 
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Fig. 8. — Epoch folded light curves derived for the 4.88 GHz (top) and 8.44 GHz (bottom) 
Stokes V radio data from TVLM 513 with higher time resolutions of 6, 7 & 8 minutes. Two 
periods of the same epoch folded light curve are shown for clarity. The large periodic peak 
in left circular polarization observed in the 4.88 GHz data has no corresponding counterpart 
in the 8.44 GHz data. In fact, the emission is mildly right circular polarized suggesting a 
reversal in polarization with frequency. 



